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Введение. Работа посвящена исследованию решения зада- 
чи трансформации форм фосфора, азота и кремния. Дан- 
ная проблема возникает при моделировании динамики 
фитопланктона в мелководных водоемах, в том числе в 
Азовском море. Модель динамики фитопланктона сфор- 
мулирована как краевая задача для системы уравнений 
диффузии-конвекции-реакции и учитывает поглощение и 
выделение питательных веществ фитопланктоном, а также 
переход питательных веществ из одного соединения в 
другое. Для расчета начальных условий и параметров 
уравнений, при которых наступает стационарный режим, 
разработано программное обеспечение, основой которого 
послужила модель, описывающая изменения концентра- 
ций фитопланктона без учета влияния течений. Данная 
модель представлена системой неоднородных обыкновен- 
ных дифференциальных уравнений. На основе разрабо- 
танного программного обеспечения экспериментальным 
образом рассчитаны начальные условия и параметры мо- 
дели динамики фитопланктона в Азовском море. 
Материалы и методы. Рассматривается трехмерная мо- 
дель динамики фитопланктона с учетом трансформации 
соединений фосфора, азота и кремния, основанная на си- 
стеме уравнений транспорта биогенных веществ. Для 
уточнения параметров модели, при которых наступают 
стационарные режимы, рассматривается случай простран- 
ственно-однородного распределения субстанций. В ре- 
зультате упрощения получена система обыкновенных 
дифференциальных уравнений, которая решена методом 
Рунге-Кутты. 

Результаты исследования. Разработано программное 
обеспечение для уточнения начальных условий и пара- 
метров модели динамики фитопланктона с учетом транс- 
формации соединений фосфора, азота и кремния. Прове- 
дены несколько численных экспериментов в предположе- 
нии, что развитие фитопланктона лимитируется един- 
ственным биогенным веществом. В результате вычисли- 
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тельного эксперимента видно, что при полученных значе- 
ниях начальных концентраций и параметров уравнений 
наступают стационарные режимы для системы обыкно- 
венных дифференциальных уравнений, описывающей 
случай пространственно-равномерного распределения 
субстанций. 

Обсуждение и заключения. В работе исследована матема- 
тическая модель трансформации форм фосфора, азота и 
кремния в задаче динамики фитопланктона. 

Получены стационарные режимы для системы обыкно- 
венных дифференциальных уравнений, для которых опре- 
делены значения параметров системы и начальные усло- 
вия. Полученные результаты могут быть использованы в 
процессе дальнейшего моделирования динамики фито- 
планктона с учетом трансформации соединений фосфора, 
азота и кремния с учетом конвекции-диффузии, солено- 
сти, температуры. 


Ключевые слова: фитопланктон, фосфор, азот, кремний, 
биоген, химико-биологический источник, уравнение кон- 
векции-диффузии-реакции, задача Коши для системы 
обыкновенных дифференциальных уравнений, стационар- 
ный режим. 
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Введение. В связи с развитием крупных городов на побережье мелководных водоемов и речных си- 


стем, впадающих в данные водоемы, участились случаи эвтрофикации. Рост водорослей в водоемах вызван 
увеличением стока азотных и фосфорных соединений с прилегающих районов суши. Каждый водный объект 
является уникальным и требует подробного исследования. Для изучения водоемов прибегают к натурным ис- 
следованиям [1], а также математическому моделированию. Не преуменьшая роли натурных экспериментов 
можно сказать, что математическое моделирование является менее затратным, а также позволяет прогнозиро- 
вать поведение экосистемы. 

Для изучения Азовского моря разработана трехмерная модель гидродинамики [2,3], включающая 
уравнения движения по трем пространственным направлениям. В работе [4] данная модель выполнена для слу- 
чая динамического перестроения геометрии расчетной области за счет приливно-отливных явлений. Исследо- 
вание точности данной модели приведено в работе [5]. В работах [6-8] приведена реконструкция экологической 
катастрофы, произошедшей в 2001 году, вызванной чрезмерной концентрацией водорослей в восточной части 
Азовского моря. В работе [9] предложены методы борьбы с заморными явлениями, возникающими в Азовском 
море. Работы [10-12] посвящены изучению динамики фито- и зоопланктона. 

Состояние вод в мелководных водоемах быстро меняется и математические модели требуют уточне- 
ния. Расчет параметров трехмерной модели динамики фито- и зоопланктона является трудоемким, поэтому 
предлагается для расчета данных параметров использовать упрощенную модель. 

Цель данной работы состоит в уточнении параметров модели динамики фитопланктона с учетом 
трансформации соединений фосфора, азота и кремния, при которых наступают стационарные режимы в пред- 
положении пространственно- однородного распределения субстанций. 

Материалы и методы. Модель основана на системе уравнений транспорта биогенных веществ 


[15, 16], вид которой для каждого модельного блока Е; имеет вид 


+и ой И. = (К отада,)+ К. , (1) 
[0] 9х д 02 ь 
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где 4, — концентрация {-ой компоненты, [мг/л]; {Е+М, М={Рь Е» Ех РО., РОР, БОР, МО, М» МН. 5%; 
{и,у, и — компоненты вектора скорости водного потока, [м/с]; А - коэффициент турбулентного обмена, [м?/с]; 
и = функция-источник биогенных веществ, [мг/л‘с]. 


В уравнении (1) индекс Г указывает на вид субстанции (таблица 1). 



































Таблица 1 
Биогенные вещества в модели динамики фитопланктона 
Номер Обозначение Название 
1 Е, зеленая водоросль СШогеНа ушеаг1$ 
2 Е› синезеленая водоросль АрВап17отепоп Ноз-адиае 
3 Е. диатомовая водоросль Зсейеопета созацит 
4 РО. фосфаты 
5 РОР взвешенный органический фосфор 
6 БОР растворенный органический фосфор 
7 МО; нитраты 
8 МО, нитриты 
9 МН. аммоний 
1 . растворенный неорганический кремний 
} ь (кремниевые кислоты) 

















Химико-биологические реакции описаны следующими уравнениями 
Е, = 9 1-К „а. —Кнь4ь —К,,Чь ‚1=1,3, 


з 
Крор — 55>К Ч, — К »Ч4рор — К рмЧрор , 
з 
Вар =- 25>К Ч + 59 ов = В Фев › 
з 
и == .5»С. те 1) Че Ааа р. > 


ти 
Я (Чо, ›Чмо, › Ч мн, 





з 3 
Кын, — 2.5» Ся (К 1) " +25» (К„ +К, ) Чл — К он, › 


и (Чо › Чмо, ›Чюн, ) Чхо, 
И (Чо › Чмо, Ч мн, ) Чхо, + Чхо, 





З 
Кхо, = 2.5»С; (К к» р Че + К оЧрн, — К Чо, › 


у. (Чо ы Чо, > Чхн, ) Чо, 
Л [бы › Чхо, ›Ч мн, ) Чко, + Чко, 





4- +К 


2зЧмо, › 


Во = 2.56, (к„„-1) 


В; = 5 Св (К -П а +5К кьЧк . 

Здесь К,„, — удельная скорость дыхания фитопланктона; К,„ — удельная скорость отмирания фито- 
планктона; К,„, — удельная скорость экскреции фитопланктона; К» — удельная скорость автолиза РОР; К»» 
— коэффициент фосфатофикации РОР; К„, — коэффициент фосфатофикации РОР; К, — удельная скорость 
окисления аммония до нитритов в процессе нитрификации; К — удельная скорость окисления нитритов до 
нитратов в процессе нитрификации, 5», 5,, 5, — нормировочные коэффициенты между содержанием №, Р, 51 в 


органическом веществе. 
Скорость роста фитопланктона определяется выражениями: 


С,, ш ет тт {. (Чл, ), № (Чл › Чло, ›Чкн, )} ? 


С; — Кук тт = (Чл ), У (Чо, › Чло, ›Чкн, р (а; )} 2 


где Кук —Щ_ максимальная удельная скорость роста фитопланктона. 
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Рис. 1. Модельная схема биогеохимической трансформации 
форм фосфора, азота и кремния 


Функции, описывающие содержания биогенов 


Чро 
— для фосфора ХЛ, (4. ) = а — | 
Чро, + РО, 
где К», — константа полунасыщения фосфатами; 
Чу 
— для кремния Х, (а ) = о. , 
Ч; Ба КУ 
где К, — константа полунасыщения кремнием; 


— для азота и. [955 › Чо, ‚Чин, )= о (Чл ›, Чко, ›Чкн, р (м. ) › 





(1) м (а + Чо, )ехр(-К „Чкн, (2) — Чун, 
м (Чо, › Чо, › ав) 2 и. (Чхн,) = > 
ба ы (Ча, ыы Чло, ) Я г Члн, 
где Ко, —_ константа полунасыщения нитратами, Кун, — константа полунасыщения аммонием, К =—= 


коэффициент ингибирования аммония. 
Для системы (1) необходимо задать векторное поле скоростей водного потока, а также началь- 


ные значения функций концентраций 4, 
а, (х,у,2,0) =а (%%2); (х,у,2) ЕС., #=0, 1ЕМ. (2) 


Будем считать, что граница » цилиндрической области С является кусочно-гладкой и 


Х=УХ, 9х, мо, где », — поверхность дна водоема, >, — невозмущенная поверхность водной среды, с 
— боковая (цилиндрическая) поверхность. Пусть и, — нормальная по отношению к » составляющая 
вектора скорости водного потока, п — вектор внешней нормали к ХУ. Для концентраций 4, будем пред- 


полагать: 


— на боковой границе: 


4, =0, на с ‚если и, <0,1=М; (3) 

94, | 

га на с, если и, 20, 1еМ; (4) 
и 


— на У, —Щ поверхность водоема: 


ам (5) 
02 


— на дне ХУ, : 
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94, 
84, ТЕ {Рь Е», Ез}, 
7 


|: =: ==, 4, = {РО., РОР, РОР, МО,, МО», МН, $1}, (6) 


где &,, =, — скорости осаждения водорослей и питательных веществ на дно. 


Стационарный режим. Рассмотрим случай пространственно-равномерного распределения суб- 
станций (фитопланктона, форм фосфора, азота и кремния), тогда каждое из уравнений (1) упрощается и 


в результате получаем следующую систему обыкновенных дифференциальных уравнений (ОДУ): 




















т =С.@-К)4, —КиьЧа — Кик» #=1,3, (7) 
ма ы УК», -К »Чы»-— КнЧ нор» 
"= Е Увек, +К „Ч ы» — К ры › 
Ч» - о С, (К -14; +К»Ч и» + К нь 
Я, =>.» С, (Ки-1) й И, ео (К + Ко 9 - Коды, 
2 _ — с, (1) Л” (№О,,М№О,,МН,) Чл И 


Л, (№О,,МО,, МН, Чло, + Чхо, 


4 
и 


Решим систему обыкновенных дифференциальных уравнений методом Рунге-Кутты [15—17]. 


[Вы -1) 4, +5,К 


врЧк 


Проведем несколько численных экспериментов, предполагая, что развитие фитопланктона зависит от 
единственного лимитирующего вещества. 
Результаты исследования. Для системы ОДУ (7) рассчитаем начальные условия и параметры 


уравнений, при которых наступают стационарные режимы. Возьмем начальные значения концентраций: 
Ч (0)= 2,5 мг/л, Ч, (0) = 2,6 мг/л, Чк (0) = 0,91 мг/л, в = 0,07 мг/л, ть (0)= 0,07 мг/л, Чро, (0) = 0,005 
мг/л, Чи, (0) = 0,11 мг/л, Чо, (0)= 0,0178 мг/л, Чо, (0) = 0,304 мг/л, аз (0) = 0,4 мг/л; коэффициенты: К, 
=2,8 сут’,К= 0,15 сут’, К„,= 0,05 сут’, К„„=0,15 сут", Крр= 0,015 сут’, Кру= 0,02 сут”, 
Крм= 0,1 сут", Кд = 0,9 сут", Кз=2,5 сут", Кри = 1,46 сут", 5р= 0,01, 5,= 0,016, 5х = 0,023, Ки, 
= 0,024, К» = 3,0, К„=2,0, К, = 3,0. 


Полученные стационарные режимы системы ОДУ (7) в предположении, что развитие фитопланк- 
тона лимитируется единственным биогенным веществом — фосфором, азотом или кремнием — изобра- 
жены на рис. 2—4 соответственно. Рис. 2 описывает влияние фосфора на развитие различных видов фи- 
топланктона, рис. 3 описывает влияние азота на развитие различных видов фитопланктона, рис. 4 опи- 


сывает влияние азота на развитие диатомовых водорослей. 
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Рис. 2. Стационарный режим системы ОДУ в предположении, что развитие фитопланктона лимитируется фосфором: 
а) зеленая водоросль (СИГ), б) синезеленая водоросль (АЁ-А), в) диатомовая водоросль (5С), г) взвешенный органический 
фосфор (РОР), д) растворенный органический фосфор (РОР), е) фосфаты (РО.) 


В результате вычислительного эксперимента видно, что при указанных выше значениях начальных 
концентраций и параметров уравнений наступают стационарные режимы для системы ОДУ (7), описывающей 
случай пространственно-равномерного распределения субстанций. Полученные значения в дальнейшем будут 
использованы для моделирования пространственно-неоднородного распределения субстанций, солености и 
температуры с учетом движения водной среды [18]. 
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Рис. 3. Стационарный режим системы ОДУ в предположении, что развитие фитопланктона лимитируется азотом: а) зеленая 
водоросль (СИГ), 6) синезеленая водоросль (АЕ-А), в) диатомовая водоросль (5С), г) аммоний (МН.), д) нитриты (МО), е) 
нитраты (№0О;) 

$1, мг/л 





Рис. 4. Стационарный режим системы ОДУ в предположении, что развитие фитопланктона (диатомовых водорослей) лими- 
тируется кремнием: а) диатомовая водоросль (5С), 6) кремний (51) 


Заключение. В данной работе исследована математическая модель трансформации форм фосфора, азо- 
та и кремния в задаче динамики фитопланктона. Рассмотрен случай пространственно-равномерного распреде- 
ления субстанций (фитопланктона, форм фосфора, азота и кремния). Система разбивается на три системы 
обыкновенных дифференциальных уравнений, каждая из которых моделирует зависимость роста фитопланкто- 
на от единственного питательного вещества. Эти системы решены методом Рунге-Кутты, получены стационар- 
ные режимы (рис. 2—4), для которых определены значения параметров системы и начальные условия. 
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Полученные результаты будут использованы для дальнейшего моделирования динамики фитопланкто- 
на с учетом трансформации соединений фосфора, азота и кремния с учетом конвекции-диффузии, солености, 
температуры. 
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